UNCLASSIFIED 

SECwR'IY  u  AsSiF  CATION  OF  ‘HIS  PAu'l 


REPORT  DOCUMENTATION  PAGE 


la  REPORT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 


2a  SECURITY  CLASSIFICATION  AUTHORITY 


2b  DECLASSIFICATION  /  DOWNGRADING  SCHEDULE 


lb  RESTRICTIVE  MARKINGS 


3  DISTRIBUTION /AVAILABILITY  OF  REPORT 

Approved  for  public  release; 
distribution  is  unlimited 


4  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 

CR  86-06 


5.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 

CR  86-06 


r.a  NAME  OF  PERFORM  NG  ORGANIZATION 

Science  Applications 
International  Corporation 


6b  OFrICE  SYMBOL 

(If  applicable) 

1-425-00-797 


7a  NA 


ME  OF  MQNITQRING  ORGANIZATION 

Naval  environmental  Prediction 


Research  Facility 


6(  A  )DRFSS  \Ufy,  State  and  ZIP  Code) 

20 o  Monteci to  Avenue 
Monterey,  CA  93940 


b  ADDRESS  (C/fy  State  and  ZIP  Code) 

Monterey,  CA  93943-5006 


83.  NAME  OF  FUNDING /SPONSORING 

ORGANIZATION 

Naval  Air  Systems  Command 


8b.  OFFICE  SYMBOL 
(If  applicable) 

AIR-330H 


9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

N00228-84-C-31 42 


8c,  ADDRESS  (City,  State  and  ZIP  Code) 

Department  of  the  Navy 
Washington,  DC  20361 


10  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM 

PROJECT 

TASK 

WORK  UNIT 

ELEMENT  NO 

NO 

NO 

ACCESSION  NO 

63207N 

7W0513 

ccoo 

DN650784 

11  TITLE  (Include  Security  Classification) 

Mid-Latitude  Wind  Probability  Program  (U) 


12  PERSONAL  AUTHOR(S) 


Jarre! 1 ,  Jerry  D. 


13a.  TYPE  OF  REPORT 

Final 


I  13b  TIME  COVERED 

FROM  _9/T2/84  to  5/1 1/86 


14  DATE  OF  REPORT  {Year,  Month,  Day)  |l  5  PAGE  COUNT 

1986  November  36 


16  SUPPLEMENTARY  NOTATION 


17 


COSATI  CODES 


FIELD 

GROUP 

SUB-GROUP 

04 

02 

18  SUBJECT  TERMS  ( Continue  on  reverse  if  necessary  and  identify  by  block  number) 

Decision  aid  Wind  probability 

Strike  probability  Mid-1 atti tudes 


19  A8STRACT  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 

A  model  to  estimate  the  probability  of  a  strike  on  a  point  by  a  mid-latitude  storm 
was  developed  and  tested.  This  model  was  based  on  a  vortex  tracking  program  operational 
for  quality  control  at  the  Fleet  Numerical  Oceanography  Center,  Monterey,  CA0  An  extension 
of  this  model  to  assess  the  probability  of  a  mid-latitude  storm  causing  30  kt  winds  also 
was  developed  and  tested.  Both  developments  followed  successful  comparable  probability 
estimation  models  for  tropical  cyclones.  Model  development  is  described  and  the  test 
methodology  and  results  are  described  and  interpreted. 


20  D'STRIBUOON 'AVAILABILITY  OF  ABSTRACT 

CHUNCLASSIFIED/UNUMITED  □  SAME  AS  RPT  □  DTIC  USERS 


22 a  NAME  OF  RESPONSIBLE  INDIVIDUAL 

Dr.  Ted  L.  Tsui 
DD  FORM  1473, 84  mar 


21  ABSTRACT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 


22b  JI LEPHONE  (Include  Area  Code)  22c  OFFICE  SYM80L 

(408)  647-4727  NEPRF  WU  6.3-14 


83  APR  edition  may  be  used  until  exhausted 
All  other  editions  are  obsolete 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


UNCLASSIFIED 


AN  (1) 
FG  (2) 
Cl  (3) 
CA  (5) 
TI  (6) 
TC  (S) 
DN  (9) 
AU  (10) 
RD  (11) 
PG  (12) 
CT  (15) 
PJ  (16) 
TN  (17) 
RN  '(18) 
RC  (20) 
DE  (,23) 


DC  ( 24 ) 
ID  (25) 

IC  (26) 
AB  (27) 


AC  (28) 
DL  ( 33 ) 
3E  (34) 
CO  (35) 


AD-A179  256 

040200 

(U) 

SCIENCE  APPLICATIONS  INTERNATIONAL  CORP  MONTEREY  CA 
Mid-Latitude  Wind  Probability  Model. 

(U) 

Final  rept.  12  Sep  84-11  May  86, 

Jarrel 1 , Jerry  D. 

Nov  1986 
38 

N00228-84-C-3142 
W0513CC 
W0513CCQ0 
NEPRF-CR-86— 06 
Unclassified  report 

*STORMS ,  * WIND ,  * WEATHER  FORECASTING,  ESTIMATES, 

LATITUDE,  MATHEMATICAL  MODELS,  MODELS,  PROBABILITY, 
QUALITY  CONTROL,  TEST  METHODS,  TRACKING,  TROPICAL 
CYCLONES,  VORTICES,  WIND  VELOCITY,  DECISION  MAKING, 
TEMPERATE  REGIONS,  NORTH  PACIFIC  OCEAN,  COMPUTER 
PROGRAMS,  CYCLONES 
(U) 

Decision  aids.  Midlatitude  storms.  Strike  probability( 

Storms),  WUDN650784 ,  PE63207N 

(U) 

A  model  to  estimate  the  probability  of  a  strike  on  a 
point  by  a  mid- latitude  storm  was  developed  and  tested. 
This  model  was  based  on  a  vortex  tracking  program 
operational  for  quality  control  at  the  Fleet  Numerical 
Oceanography  Center,  Monterey,  CA.  An  extension  of  this 
model  to  assess  the  probability  of  a  mid-latitude  storm 
causing  30  kt  winds  also  was  developed  and  tested.  Both 
developments  followed  successful  comparable  probability 
estimation  models  for  tropical  cyclones.  Model 
development  is  described  and  the  test  methodology  and 
results  are  described  and  interpreted.  Keywords: 
Decision  aid,  Strike  probability,  Wind  probability, 
Mid-latitudes,  Wind,  Wind  velocity. 

(U) 

01 

F 

417277  _ _ 


CONTENTS 


r 

1.  Introduction  .  1 

2.  Strike  Probability  .  3 

3.  Wind  Probability  Model  .  8 

4.  Independent  Testing  .  14 

4.1  Verification  Data . 14 

4.2  Test  Results  . . 15 

5.  Summary  and  Recommendations  .  20 

References  .  21 

Appendix  A  --  User's  Manual  for  Program  WNPPSG  .  .  A-l 

Distribution  .  Dist-1 


s* 


i 


NAVENVPREDRSCHFAC  CR  86-06 


I 


t 


U  Wax  ' 

S*plv  r-mm 

MONTERtV,  CAUFORBIfl 


NAVENVPREDRSCHFAC 
CONTRACTOR  REPORT 
CR-86-06. 


MID-LATITUDE  WIND  PROBABILITY  MODEL, 


Jerry  D,  Jarrell 

Science  Applications  International  Corp. 
Monterey,  CA  93940 


Contract  No.  N00228-84-C-3142 


NOVEMBER  1986 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  IS  UNLIMITED 


Prepared  For: 

NAVAL  ENVIRONMENTAL  PREDICTION  RESEARCH  FACILITY  , 
11  MONTEREY,  CALIFORNIA  93043-5006 


« 


1.  Introduction 


The  concept  of  using  probabilities  to  extend  the  useful 
life  of  storm  forecasts  was  attempted  nearly  thirty  years  ago 
by  Kimball  (1958)  and  later  Appleman  (1962).  Both  used 
graphical  means  to  estimate  the  probability  of  typhoon  winds 
around  U.S.  Air  Bases  in  the  Far  East.  Practical  use  of 
probabilities  began  when  Jarrell  (1978)  developed  computeriz¬ 
ed  methods  to  sum  the  tropical  cyclone  probabilities  over 
time.  The  best  known  of  these  programs  is  the  strike  proba¬ 
bility  ( STRIKP )  which  the  Navy,  Air  Force  and  National 
Weather  Service  use  to  assess  the  likelihood  of  a  tropical 
cyclone  passing  close  enough  to  a  base  or  city  to  cause  major 
damage.  Therein  lies  the  definition  of  strike  —  that  the 
cyclone  passed  close  enough  to  cause  major  damage.  STRIKP  is 
driven  by  point  forecasts  issued  by  a  tropical  cyclone  warn¬ 
ing  center. 

The  idea  of  adapting  the  STRIKP  model  to  provide  comp¬ 
arable  sophistication  to  the  mid  latitude  forecasters  has 
been  hampered  by  the  lack  of  centralized  consistent  point 
forecasts  on  mid  latitude  cyclones.  However,  point  forecasts 
in  a  useable  format  became  available  with  the  development  of 
the  Quality  Control  Vortex  Tracking  Program  (QCVTP)  (see  Tsui 
et  al.  1982) . 

The  purpose  of  this  study  is  to  demonstrate  that  the 
•  proven  strike  probability  concept  can  be  successfully  operat¬ 

ed  using  the  QCVTP  data  bases.  It  should  be  noted  that  a 
x  statistical  proof  of  the  concept  is  not  intended. 
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The  concept  of  a  mid-latitude  wind  probability  model  was 
formulated  by  SAIC  under  NEPRF  contract  No.  N00228-84-C-3142 . 
The  need  for  a  mid-latitude  wind  probability  model  is  to  aid 
ship  captains  and  other  decision  makers  in  evading  damaging 
winds  from  mid-latitude  or  extratropical  cyclones.  The 
cyclones  are  forecast  but  each  forecast  element  (e.g.,  track, 
intensity  and  size)  is  subject  to  error.  Wind  probability 
helps  allow  for  error  by  preventing  overstimulation  at  points 
along  the  forecast  track  and  understimulation  at  points  off 
the  track  which  are  also  significantly  threatened. 

The  model  accepts  QCVTP  type  storm  input  data  (low  pres¬ 
sure  center  latitude,  longitude,  and  pressure  deficit,  as 
well  as  the  radius,  ellipticity  and  angle  of  the  major  axis 
of  the  zero  anomaly  countour  (ZAC))  at  analysis  time  and  up 
to  5  prognosis  times  (12-60  hrs).  The  output  consists  of 
strike  and  gale  force  wind  probabilities  in  both  instantane¬ 
ous  and  elapsed  time  modes. 


r 
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2. 


Strike  Probability 


y 


* 


The  tropical  cyclone  working  definition  of  a  "strike"  is 
that  the  cyclone  passes  within  an  area  defined  by  semicircles 
of  75  and  50  n  mi  with  the  larger  being  on  the  left  in  the 
northern  hemisphere  (facing  direction  of  motion).  This 
allows  for  a  greater  extent  of  strong  winds  on  the  right  side 
(see  Jarrell  and  Brand,  1983).  This  definition  appears 
inappropriate  here  since  there  is  a  large  area  of  fairly 
light  winds  near  the  center  of  an  extratropical  cyclone  and 
then  the  strong  winds  occur  in  a  broad  region  often  extending 
hundreds  of  miles  from  the  center. 


Figure  1.  Scatter 
diagram  of  Oct-Dec 
1984  North  Pacific 
extratropical  cy¬ 
clones  minimum  sea 
level  pressure  vs 
pressure  at  the 
zero  anomaly  con¬ 
tour  (ZAC).  The 
line  indicates 
the  demarcation 
between  strike 
area  definitions 
using  the  1000 
isobar  below  the 
line  and  one  half 
ZAC  above  the 
line. 
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It  was  first  thought  that  perhaps  a  more  reasonable 
definition  would  be  entry  into  the  area  of  depressed 
pressure,  or  within  the  zero  anomaly  contour  (Tsui  et  al . 
1982).  This  area  is  quite  large— often  5  grid  lengths  and 
greater  in  radius  (a  grid  length  is  .206  n  mi  at  60°N); 
consequently,  being  struck  would  then  mean  little  more  than 
being  in  an  area  of  disturbed  weather.  The  strike  definition 
that  was  adopted  is  the  storm  passing  close  enough  so  the 
station  falls  within  the  1000  mb  isobar.  For  storms  with  a 
high  central  pressure  so  that  the  area  within  the  1000  mb 
isobar  is  small  or  nonexistent,  this  strike  area  is  limited 
on  the  low  side  by  half  the  radius  of  the  zero  anomaly 
contour  (measured  in  grid  lengths).  The  demarcation  between 
the  two  definitions  of  the  strike  area  is  shown  as  a  straight 
line  in  figure  1.  Strike  probability  may  therefore  be 
defined  as: 


ps  =  p(d  <  Rs> 


where  Rs  is  the  strike  radius;  RS=KRZ ,  and  Rz=the  ZAC 


radius.  The  multiplier  K  is  defined  as  the  maximum  of 
[(central  pressure  -  1000 )/( pressure  deficit)  or  0.5]  and  d 
is  the  passing  distance.  Now  since  both  d  and  Rs  are 
random  variables  (although  both  are  forecast,  there  is 
uncertainty  in  their  actual  values)  it  is  necessary  to 
restate  the  above  expression  as: 


1 


Ps  =  |  P ( R 


s  2. 


d  =  r )  dr  =  J  P(RS  >  r|d  =  r)  P(d  =  r)dr. 


In  this  initial  model  the  .assumption  is  made  that  size  errors 
are  independent  of  position  errors  and  hence  that 
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P  ( Rs  >  r|  d  =  r)  =  P(RS  >  r). 

Thus  the  expression  to  be  integrated  becomes: 

Pg  =  ^fp(Rs  >  r)  P  (  d  =  r )  dr 

The  expression  P(Rs>r)  was  evaluated  by  creating  a 
probability  contingency  table  of  the  radius  of  the  ZAC 
observed  (in  the  QCVTP  data)  versus  that  forecast  (see  Table 
1).  The  probability  term  evaluated  by  the  use  of  a 
contingency  table  is  being  estimated  on  the  basis  of 
climatological  experience.  This  could  have  been  accomplished 
by  parameterizing  the  distribution;  however,  from  the 
standpoints  of  sample  size  and  convenience  of  computation  the 
contingency  table  is  preferable.  Thus  the  first  expression 
P(Rs>r)  can  be  viewed  as  P(RS  —  r|  Rf)  P(Rf)  where 

P(  Rf )  is  the  probability  of  forecast  Rf  being  made; 
P(Rf)  =  1.0,  since  the  forecast  is  known.  Table  values  are 
given  for  r  =  <2 . 5 , 2 . 5-3 . 0  ,  3. 0-3. 5  ...  >  5.5  in  multiples 
of  grid  lengths.  Since  r  takes  on  these  discrete  values, 
the  above  integral  is  treated  as  a  summation  and  the  term 
P(d  =  r)  is  evaluated  as  the 

P(r  -  .25  <  d  _<  r  +  .25) 

which  is  a  ring  about  the  point  of  interest  of  width  one  half 
grid  length.  This  term  is  position  dependent.  Required  is 
the  probability  that  the  storm  center  falls  within  concentric 
rings  about  the  point  of  interest .  If  the  rings  were  center¬ 
ed  on  the  storm  forecast,  then  a  contingency  table  approach 
could  again  be  used  at  considerable  savings  in  program  com¬ 
plexity  and  execution  time.  The  most  reasonable  way  to  eval¬ 
uate  this  term  is  by  fitting  a  parametric  distribution  to 
historical  forecast  errors.  The  distribution  is  then  inte- 
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grated  over  rings  with  a  variable  center  location  relative  to 
the  distribution  center.  A  bivariate  normal  distribution  has 
typically  been  used  to  fit  such  errors  because  it  has  charac¬ 
teristics  which  match  those  of  the  error  population  (posi¬ 
tively  and  negatively  unbounded  with  a  central  tendency  in 
two  components).  See  Crutcher  et  al^  1982  for  a  discussion 
of  this  distribution  and  forecast  errors.  The  goodness  of 
fit  has  not  been  determined  in  this  case. 


Table  1.  The  probability  of  a  particular  mean  ZAC  radius 
being  exceeded  in  the  verifying  analysis  given  the  twenty- 
four  hour  forecast.  Similar  tables  are  used  for  other  fore¬ 
cast  intervals.  The  column  under  the  0.0  grid  lengths  head¬ 
ing  would  be  100%  except  for  dying  storms  whose  frequency  is 
shown  under  the  column  NONE. 


24-HOUR  TABLE  OF  EXCEEDANCE  PROBABILITIES  OF  R 


0.0 

2.5 

3.0 

3.5 

4.0 

4.5 

5.0 

5.5 

NONE 

SAMPLE 

FORECAST 

<2.5 

72. 

19. 

10. 

7. 

3. 

0. 

0. 

0. 

28. 

89 

2. 5-3.0 

78. 

51. 

22. 

7. 

3. 

0. 

0. 

0. 

22. 

68 

3. 0-3. 5 

90. 

80. 

59. 

26. 

13. 

3. 

1. 

0. 

10. 

93 

3 .5-4.0 

89  . 

85. 

70. 

51. 

22. 

6. 

3. 

1. 

11. 

94 

4.0-4 .5 

86. 

85. 

77. 

68. 

41. 

14. 

5. 

1. 

14. 

73 

4. 5-5.0 

85. 

82. 

78. 

72. 

51. 

29  . 

12. 

3. 

15. 

68 

5. 0-5. 5 

88. 

84. 

77. 

74. 

67. 

53. 

28. 

12. 

12. 

43 

>5.5 

89. 

87  . 

87. 

80. 

71. 

64. 

50. 

36. 

11. 

70 

Time  integration 

is  handled 

just 

as 

it  is 

in 

the 

tropical 

programs  except  for  the  following: 


a.  Given  the  tropical  input  from  warnings,  all  fore¬ 
cast  sequences  run  from  tau  =  0  to  some  last  tau, 
usually  72  hours.  With  mid-latitude  cyclones,  the 
low  frequently  appears  first  on  a  prog  chart. 

b.  Points  are  frequently  threatened  by  more  than  one 
low  center  over  the  60  hour  period  of  integration. 
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Figure  2b.  Elapsed  time  strike  probability  (units  are  tens 
of  % )  . 

Figure  2.  24-hour  probability  forecast  for  1200  GMT  8 
November  1984,  North  Pacific  Ocean. 


Probabilities  from  different  storms  are  summed  as 
if  independent. 

P ( A  or  B)  =  P ( A )  +  P ( B )  -  P ( A  and  B) 
and  P( A  and  B)  =  P(A)P(B) 

c.  All  time  integrated  probabilities  are  given  for  the 
entire  60  hour  period.  Exception  (a)  (above) 
causes  probabilities  to  suddenly  appear  on  a  time 
line  whereas  exception  (b)  causes  humps  on  the  time 
line  unlike  the  smooth  increase  to  a  peak  then  a 
decrease  over  time  as  seen  in  the  tropical 
probability  curve. 

Sample  strike  probabilities  are  shown  in  figure  2. 


Figure  2a.  Strike  probability  (units  are  tens  of  %) 
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3. 


Wind  Probability  Model 


The  wind  probability  model  is  designed  to  produce  the 
probability  of  gale  force  winds  ()>28  kt)  although  provision 
was  retained  to  include  one  other  wind  speed  (e.g.,  storm 
force  or  48  kt )  .  Both  instantaneous  and  elapsed  time  wind 
probabilities  will  be  provided  in  the  same  format  as  strike 
probabilities . 

The  form  of  the  wind  probability  is 
Jp(VD  >  28|  x , y )  P (  x , y )  dxdy 

where  VQ  =  the  wind  speed  at  point  x0,y0,  the  point  of 
interest;  x,y  is  the  location  of  the  storm  center.  Clearly 
P(VG  >  28|  x,y)  is  highly  dependent  on  the  location  of  point 
x0,y0  relative  to  point  x,y.  What  is  needed  is  a  model 
to  relate  the  wind  probability  to  locations.  The  relation¬ 
ship  of  wind  to  location  relative  to  the  storm  center  is  a 
function  of  storm  size,  pressure  depth,  the  shape  of  the 
pressure  profile  and  the  assymetry  of  the  storm. 

In  an  analysis  mode  these  parameters  can  be  measured  and 
converted  to  good  wind  estimates  by  use  of  the  gradient  or 
geostrophic  wind  equations.  This  also  applies  to  the  prog¬ 
nostic  mode  except  that  the  various  parameters  are  of  unknown 
accuracy;  hence  the  forecast  winds  are  of  unknown  value.  In 
a  probabilistic  mode  one  wants  to  find  the  sets  of  parameters 
which  produce  winds  of  interest,  for  example,  over  28  kt  or 
gale  force  winds.  The  joint  probability  of  these  sets  occur¬ 
ring  would  then  be  estimated.  The  more  variables  (e.g., 
size,  minimum  pressure,  profile  shape,  as  well  as  direction 
and  distance  of  the  point  of  interest  from  the  storm) 
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involved  in  each  alogrithm  the  more  complex  the  problem. 
Also  additional  variables  cause  additional  stratification  of 
the  available  data  base;  thus  the  size  of  the  data  base  in 
effect  limits  the  complexity  of  the  joint  probability 
solution . 

In  the  tropical  models,  the  wind  field  was  assumed  to  be 
totally  defined  by  three  parameters,  maximum  wind,  a  scaling 
radius  and  forward  speed.  Forward  speed  had  a  rather  minor 
role  and  was  assumed  to  be  perfectly  forecast.  This  left 
only  two  parameters  to  describe  the  wind  field.  Their  joint 
probability,  given  the  forecast,  was  readily  determined  from 
climatological  forecast  data. 

There  are  no  similar  models  for  extratropical  cyclones. 
It  is  assumed  here  that  the  assymetric  radial  profile  of 
winds  is  a  function  of  the  maximum  pressure  gradient  and 
orientation  of  the  major  ellipse  axis.  The  maximum  pressure 
gradient  (in  a  mean  sense)  can  be  computed  by  the  pressure 
differential  divided  by  the  length  of  the  minor  axis. 
Consequently  storms  were  classed  by  maximum  gradient  vs  angle 
of  orientation  of  the  major  axis. 

The  distribution  of  these  classes  is  shown  in  the 
scatter  diagram  of  Figure  3.  Presumably  storms  with  the 
strongest  mean  gradient  have  the  strongest  winds.  Next 
contingency  tables  of  forecast  classes  versus  observed 
classes  were  computed;  thus  the  probability  of  a  particular 
class  occurring,  given  a  forecast,  could  now  be  estimated 
(see  Table  2).  The  spatial  distribution  of  winds  was  now 
derived  separately  for  each  storm  class. 
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Table  2.  The  probability  of  a  storm  becoming  a  member  of 
each  size-depth  class  after  24  hours  given  its  particular 
forecast  class.  Percentages  left  to  right  in  a  row  would  add 
4  to  100%  except  for  dying  storms. 


24-HOUR 

CLASS  CONTINGENCY 

OBSERVED  CLASS 

TABLE 

SAMPLE 

FCST 

1 

2  3 

4  5  6  7 

8  9 

SIZE 

CLASS 


1 

22. 

19  . 

6. 

10. 

8 

4. 

11. 

5 

1. 

79 

2 

9  . 

40 . 

13  . 

6. 

11 

4. 

4. 

4 

4. 

53 

3 

3. 

26. 

37. 

3. 

3 

0  . 

5. 

3 

5. 

38 

4 

7. 

7. 

4. 

26. 

14 

4. 

11. 

8 

1. 

114 

5 

2. 

5  . 

3. 

11. 

32 

.  18. 

5. 

1 

2. 

87 

6 

3. 

2. 

8. 

12. 

10 

.  31. 

0. 

2 

5. 

59 

7 

9  . 

3. 

2. 

12. 

8 

0. 

26. 

20 

3  . 

65 

8 

3. 

12  . 

1. 

7. 

12 

6. 

7. 

19 

.  15 . 

68 

9 

2. 

7. 

9  . 

2. 

4 

7. 

9. 

22 

.  24 . 

46 

The  approach  was 

to 

compute  winds  on 

a 

grid 

centered  on 

the 

cyclone 

• 

This 

is 

a 

special 

symmetrical 

grid  with 

concentric 

circles 

of 

grid 

points 

increasing 

in  number 

outward  so 

that 

all 

grid  spaces  contain  equal 

areas  ( see 

Figure  4).  Stored  at  grid  points  for  each  class  storm  is  the 
probability  that  gale  force  winds  would  be  observed  at  the 


center  or  grid  space  #1  if  a  storm  (of  the  given  class)  were 
located  in  the  grid  space  in  question  (see  Figure  5). 
Assymmetry  is  considered  at  the  expense  of  a  somewhat  complex 
model  since  there  are  now  71  grid  points  for  each  integration 
step.  However,  because  the  grid  spaces  are  equi-area,  the 
computation  of  the  probability  of  the  storm  occupying  a  grid 
space  is  simplified. 


* 
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Figure  4.  Equiarea  grid  which  was 
placed  over  moving  cyclone  to 
extract  wind  speeds  for  evaluation 
of  conditional  wind  probabilities. 
The  radius  of  the  central  circle 
is  about  85  n  mi  and  the  outer 
circle  720  n  mi.  Numbers  are  grid 
identifying  numbers. 


Figure  5.  Grid  of  gale  force  wind 
probabilities  for  class  three 
N-NE  and  strong)  cyclones  given 
in  whole  percent.  This  is  based 
only  on  October  1984  -  April  1985 
data . 


With  this  grid  oriented  model  the  operative  wind  prob¬ 
ability  expression  becomes: 


P(V>Vg»c-K,  g-lj  c=k|  g=L)  P(g=I^  forecast 
x,y)  P(c=k|  forecast) 


where  the  first  product  term  is  the  probability  that  the  wind 
equals  or  exceeds  gale  force  given  that  a  storm  of  class  =  K 
is  located  in  grid  space  L.  The  second  product  term  is  the 
probability  that  a  storm  occupies  grid  space  L  given  the 
forecast  for  this  time  x,y.  The  last  term  is  the  probability 
that  the  storm  is  of  class  K  given  the  forecast  class.  An 
assumption  of  independence  between  the  last  two  terms  is 
implicit . 
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The  resulting  probabilities  are  reasonable  in  appear¬ 
ance.  Sample  24-hour  wind  probabilities  are  shown  in  figure 
6. 


Figure  6a. 
%)  . 


Gale  force  wind  probability  (units  are  tens  of 


IreUtJ„sbof  %E,1.apSed  time  force  wind  probability  (units 


Figure  6.  24-hour  wind  probability  forecast 
November  1984,  North  Pacific  Ocean. 


for  1200  GMT  8 
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4. 


Independent  Testing 


An  independent  data  set  of  FNOC  analysis  and  prognosis 
surface  pressure  fields  was  obtained  for  December  1983 
through  March  1984.  Forecast  intervals  were  from  12  to  48 
hours  by  12  hour  increments.  The  QCVTP  program  was  run  on 
these  fields  for  the  North  Pacific  Ocean,  providing  the  basis 
for  the  strike  and  wind  probability  computation.  To  maintain 
approximate  serial  independence  between  fields,  every  third 
12  hourly  date-time-group  was  selected  for  probability  compu¬ 
tations.  In  all,  56  date-time  groups  were  chosen  which 
involved  99  different  cyclones  which  appeared  in  an  average 
of  2.2  prognosis  series  each  (there  are  244  12  hourly  fields 
in  the  December-March  period:  however,  several  fields  were 
missing )  . 

4.1  Verification  Data 


Since  probabilities  would  be  computed  on  prognosis 
series  of  48  hours  long  initiated  every  36  hours,  all  analy¬ 
ses  would  be  used  as  verification  data  (except  for  the  first 
few).  For  the  purpose  of  strike  verification,  each  cyclone 
was  examined  against  each  grid  point  over  the  Pacific  sector 
of  the  FNOC  63x63  northern  hemisphere  grid.  A  strike  was 
considered  to  have  occurred  if  the  grid  point  fell  within  an 
elliptical  strike  area  surrounding  the  cyclone.  This  strike 
area  was  defined  as  an  ellipse  approximating  the  lOOOmb  con¬ 
tour  or  within  1/2  the  distance  out  to  the  zero  anomaly  con¬ 
tour,  which  ever  was  greater.  A  file  was  created  showing 
which  grid  points  had  been  struck  for  each  available  analysis 
date-time-group.  Gale  force  winds  were  verified  in  much  the 
same  way  except  that  geostrophic  winds  were  calculated  at 
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grid  points.  When  80%  of  the  geostrophic  wind  exceeded  gale 
force  (28kt)  it  was  considered  that  gale  force  winds  had  been 
observed  at  the  grid  point.  The  factor  of  80%  is  used  to 
compensate  for  the  fact  that  geostrophic  winds  generally  over 
estimate  actual  surface  winds  and  particularly  so  around 
cyclones.  As  with  strikes,  a  file  was  created  showing  which 
grid  points  had  received  gale  force  winds. 

With  the  strike  and  wind  verification  data  in  place, 
each  selected  date-time  was  processed  by  computing  strike  and 
wind  probabilities  for  each  grid  point.  These  were  then 
categorized  into  5%  cells  and  a  running  sum  of  the 
probabilities,  the  verification  (0  if  not  struck,  1  if  struck 
and  similiar  for  winds)  and  case  count  were  maintained  for 
each  cell.  Comparisons  were  made  of  the  expected  occurrence 
of  strikes  or  gale  force  winds  versus  the  observed 
occurrences . 

4.2  Test  Results 


Figures  7  through  10  show  the  results  of  these 
comparisons.  The  straight  line  from  lower  left  to  upper 
right  represents  the  expected  frequency  of  occurrence.  The 
asterisk  represents  the  actual  frequency  of  occurrence.  The 
vertical  I-beam  represents  99%  confidence  limits  around  the 
expected  frequencies  based  on  a  "t"  test.  Visually  the 
agreement  between  observed  strike  frequencies  and  expected 
strike  frequencies  is  quite  good.  However,  it  should  be 
noted  that  far  too  many  (much  more  than  1%)  of  the  observed 
values  fall  outside  the  99%  confidence  limits.  The  visual 
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Figure  7.  Comparison  of  expected  (diagonal  line)  vs  the 
actual  (*)  strike  frequencies.  I-beams  represent  99%  confi¬ 
dence  limits  (t-test)  about  the  expected  frequency.  Both  the 
ordinate  and  abscissa  are  frequencies  (0  to  1). 


fit  confirms  that  the  strike  probability  is  a  good  estimate 
of  the  actual  strike  frequency,  while  the  excessive  number  of 
cases  falling  outside  the  99%  confidence  limits  points  out 
that  there  are  limitations  on  the  applicability  of  at  least 
some  of  the  model  assumptions.  For  example,  in  the  instan¬ 
taneous  mode  there  is  apparently  a  tendency  to  underestimate 
strike  probabilitites .  This  can  be  attributed  to  better  FNOC 
model  forecasting  in  the  independent  cases  than  in  the  devel¬ 
opment  sample.  It  would  not  be  surprising  that  such  changes 
would  occur  either  because  the  models  change  (note  the  inde¬ 
pendent  set  predates  by  one  year  the  development  set)  or 
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Figure  8.  Comparison  of  expected  (diagonal  line)  vs  the 
actual  (*)  elapsed  time  strike  frequencies.  I-beams  repre¬ 
sent  99%  confidence  limits  (t-test)  about  the  expected  fre¬ 
quency.  Both  the  ordinate  and  abscissa  are  frequencies  (0  to 


simply  due  to  a  natural  difference  in  the  two  rather  limited 
data  sets.  Notice  that  in  the  elapsed-time-strike-probabili- 
ties  there  is  a  tendency  to  significantly  overestimate  proba¬ 
bilities  which  runs  counter  to  the  previously  mentioned  prob¬ 
lem  and  consequently  the  agreement  between  expected  and 
observed  improves  over  time.  Overestimation  in  the  summation 
process  is  caused  by  a  faulty  parameterization  of  the  inde¬ 
pendence  of  strikes  in  consecutive  time  steps.  This  would  be 
an  overestimate  of  independence  perhaps  because  large  time 
steps  allow  substantial  storm  movement  between  probability 
computations . 

The  comparison  for  gale  force  winds  is  a  different  story  (see 
figures  9  and  10).  Visually  it  is  apparent  that  there  is 
little  predictive  skill  in  the  wind  probabilities  except  at 
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0  36HR.  WINDP  1 
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0  24HR.  WINDP  1 


Figure  9.  Comparison  of  expected  (diagonal  line)  vs  the 
actu3l  (*)  gale  force  wind  frequencies.  I— beams  represent 
99%  confidence  limits  (t-test)  about  the  expected  frequency. 
Both  the  ordinate  and  abscissa  are  frequencies  (0  to  1). 


low  probability  levels  and  even  then  they  are  biased  on  the 
low  side.  It  is  not  likely  that  this  poor  agreement  can  be 
attributed  to  one  cause,  but  the  lack  of  a  good  model  of  the 
spatial  distribution  of  winds  around  extratropical  cyclones 
is  doubtless  a  major  contributing  factor.  QCVTP  by  its  meth¬ 
odology  of  fitting  ellipses  to  low  pressure  aras  is  not  a 
good  estimator  of  wind  fields.  The  low  side  bias  may  be 
attributable  to  using  80%  of  the  geostrophic  wind  for  verifi¬ 
cation  and  indeed  using  75%  of  that  wind  essentially  elimates 
that  bias,  but  of  course  the  agreement  with  the  high  proba¬ 
bilities  in  comparably  aggravated  since  they  are  generally 
over  forecast.  One  of  the  problems  which  plagues  the  wind 
probabilities  is  that  there  tends  to  be  a  very  large  central 


18 


4 


Figure  10.  Comparison  of  expected  (diagonal  line)  vs  the 
actual  (*)  elapsed  time  gale  force  wind  frequencies.  I-beams 
represent  99%  confidence  limits  (t-test)  about  the  expected 
frequency.  Both  the  ordinate  and  abscissa  are  frequencies  (0 
to  1 )  . 


area  within  an  extratropical  cyclone  which  is  characterized 
by  light  winds.  For  this  reason  very  high  strike  probabili¬ 
ties  are  associated  with  light  winds.  These  also  tend  to  be 
associated  with  high  wind  probabilities  as  can  be  seen  in 
figures  9  and  10. 
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5. 


Summary  and  Recommendations 


A  realistic  model  has  been  developed  and  tested  for  mid¬ 
latitude  strike  probabilities.  A  comparable  model  for  wind 
probabilities  proved  to  be  much  less  reliable  and  not  usable. 
The  definition  of  strike  does  not  necessarily  imply  severe 
weather  from  the  standpoint  of  naval  operations,  but  it  does 
imply  that  one  of  the  necessary  conditions  required  for 
severe  weather,  the  presence  of  a  cyclonic  storm,  has  been 
met.  With  subjective  screening  of  the  severity  of  a  storm, 
STRIKP  can  be  used  as  an  early  warning  precursor  of  severe 
weather  --  perhaps  the  basis  of  a  storm  or  gale  condition  III 
or  IV  (gales  are  possible  within  48  or  72  hours  respective¬ 
ly).  Similarly,  for  ship  routing,  a  high  strike  probability 
may  not  be  suitable  justification  for  definitive  evasive 
actions  but  rather  to  position  oneself  favorably  for  evasion 
should  it  be  warranted  later. 
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APPENDIX  A 

user's  MANUAL  FOR  PROGRAM  WNPPSG 
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1.0  Introduction 

The  Pacific  mid-latitude  strike  and  wind  probablilty  model  was 
developed  by  5AIC  under  NEPRF  contract  No.  N00228-84-C-3142.  The 
reason  for  developing  a  mid-latitude  cyclone  probability  model  is  to  aid 
ship's  captains  and  other  decision  makers  in  evading  dangerous  conditions 
in  mid-latitude  or  extratropical  cyclones  Track,  intensity,  and  size  of 
these  cyclones  are  forecast,  but  each  forecast  element  is  subject  to  error. 
Probablity  helps  allow  for  error  by  preventing  overstimulation  at  points 
along  the  forecast  track  and  understimulation  at  points  off  the  track 
which  are  also  significantly  threatened. 

The  model  accepts  5EIS  output  storm  data  ,(low  pressure  center 
latitude,  longitude,  pressure  deficit,  as  well  as  the  radius,  ellipticity  and 
angle  of  the  major  axis  of  the  zero  anamoly  countour  (ZAC))  at  analysis 
time  and  up  to  5  prognosis  times  (12-60  hrs),  as  input.  The  output 
consists  of  strike  and  gale  force  wind  probabilities  in  both  instantaneous 
and  time  elapsed  modes. 
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2.0  Program  Operation 


The  Model  software  is  called  WNPPSG  (^iNd  Probability,  Polar 
Sterographic  Grid).  WNPPSG  accepts  real  time  output  from  the  5EIS 
program  as  input.  The  files  from  SEIS  are:  ANSTRM,  FCSTRM,  and  QCDATE. 
QCDATE  is  the  current  date  on  which  SEIS  was  run.  ANSTRM  is  in  Update 
format  and  is  updated  by  the  JCL..  ANSTRM  consists  of  low  pressure 
center  latitude,  longitude,  and  pressure  deficit,  as  well  as  the  radius, 
ellipticity  and  angle  of  the  major  axis  of  the  zero  anamoly  countour  (ZAC) 
at  analysis  time  and  up  to  5  prognosis  times  (12-60  hrs).  FCSTRM 
consists  of  the  same  parameters  of  storms  that  are  forecast  to  occur  but 
which  have  not  yet  been  observed. 


WNPPSG  also  reads  files  PDLTR  and  CLSTAT.  PDLTR  is  a  file  of 
bivariate  normal  probabilities  and  CLSTAT  is  a  file  of  climatological 
statistics  on  the  frequency  of  occurrence  of  storms  of  various  radii  and  of 
storm  classes  based  on  orientation  versus  maximum  pressure  gradient. 

Output  is  a  series  of  plots  of  strike  and  gale  force  wind  probabilities 
in  both  instantaneous  and  time  elapsed  modes. 

The  Technical  aspects  of  this  program  are  described  in  the  final 
report  for  Contract  N00228-84-C-3 1 42. 
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The  following  are  some  examples  of  the  real  time  input  files  used  by 
WNPPSG. 
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Figure  1.  A  segment  of  the  File  ANSTRM,  16  storms  have  been  omitted. 
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File  FCSTRM: 
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Figure  2.  A  segment  of  the  File  FCSTRM,  2  storms  have  been  omitted. 
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JCL  for  the  Program  WNPPSG. 


xxxxl,ST5PC,5J,P3,EC  1 50.  WNPPSG 
ATT  ACH.QCPAC,  I  D=XJ. 
ULIB(QCPAC,ANSTRM,FC5TRM,QCDATE) 

RETURN,QCPAC. 

UPDATE,F,D,L-0,C-0,I-EMPTY,P-ANSTRM,S-TAPE3. 

REWIND,FCSTRM. 

C0PY,FCSTRM,TAPE4 
ATT  ACH,CLST  AT,  I  D”JL. 

REWIND,CLSTAT. 

COPY,CLST  AT,T  APE5. 

RETURN, CLSTAT. 

ATTACH,  PDLTR,  I  D=JL. 

REWIND, PDLTR. 

COPY, PDLTR, TAPE6. 

RETURN, PDLTR. 

REWIND,QCDATE. 

COPY,QCDATE,T  APE8. 

RETURN,  AN5TRM,FCSTFtf1,QCD  ATE. 

L I BR  ARY,*FNWCL  I B. 

MAP, ON. 

FTN,R-2. 

LGO. 

MAP,OFF. 

LIBRARY(*FNWCOVL) 

ROUTE,  PLOT,  DC-PV,FID-xxxx,TID«C,BIN-xx,DEF,CO-xxxxxx. 
0NSW,6. 

VARITNY. 

*EOR 
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Figure  3.  24  Hour  Elapsed  Time  Gale  Force  Wind  Probability  for  the  Pacific. 
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Figure  4.  24  Hour  Gale  Force  Wind  Probability  for  the  Pacific. 
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Figure  5.  24  Hour  Elapsed  Time  Strike  Probability  for  the  Pacific. 
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Figure  6.  24  Hour  Strike  Probability  for  the  Pacific. 
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